Decoupling of the S — 1/2 antiferromagnetic zig-zag ladder with anisotropy 
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The spin- 1/2 antiferromagnetic zig-zag ladder is studied by exact diagonalization of small systems 
in the regime of weak inter-chain coupling. A gapless phase with quasi long-range spiral correlations 
has been predicted to occur in this regime if easy-plane {XY) anisotropy is present. We find in 
general that the finite zig-zag ladder shows three phases: a gapless collinear phase, a dimer phase 
and a spiral phase. We study the level crossings of the spectrum, the dimer correlation function, the 
structure factor and the spin stiffness within these phases, as well as at the transition points. As 
the inter-chain coupling decreases we observe a transition in the anisotropic XY case from a phase 
with a gap to a gapless phase that is best described by two decoupled antiferromagnetic chains. 
The isotropic and the anisotropic XY cases are found to be qualitatively the same, however, in the 
regime of weak inter-chain coupling for the small systems studied here. We attribute this to a finite- 
size effect in the isotropic zig-zag case that results from exponentially diverging antiferromagnetic 
correlations in the weak-coupling limit. 
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X 



I. INTRODUCTION 



Antiferromagnetic ladder systems have attracted much 
interest recently.El On the theoretical side they interpct. 
late between the well studied antiferromagnetic chairo 
and two dimensional antiferromagnets.u The evolution 
between ID and 2D spin- 1/2 antiferromagnetism is not 
necessarily smooth, however. In particular, the n-leg lad- 
der shows a remarkable alternating property in the spec- 
trum as the number of legs is even or odd.u The spectrum 
has a gap for an even number of legs while it is gapless for 
an odd number of legs. This is similar to the difference 
between integer (spin-gap) and half-odd-integer (gapless) 
spin chainsEl. In the limit of strong-coupling between the 
two chains the two-leg ladder is essentially composed of 
weakly interacting singlets that form across the rungs. 
The lowest excitation is the promotion of a rung singlet 
to a triplet with an excitation energy of the order of the 
inter-chain coupling. This spin gap remains nonvanishing 
even for small inter-chain coupling due to, the fact that a 
single antiferromagnetic chain is critical. □ In the case of 
purely Ising coupled chains the gap appears for all values 
of nil 

The antiferromagnetic zigzag ladder has also attracted 
interest recently, particularly in the context of exper- 
imental systems with low-dimensional magnetic struc- 



tures like that of C^CuCZja. It is also interesting from 
a theoretical point of view because it is a frustrated sys- 
tem (see Fig. 1). Indeed, the zig-zag ladder is equivalent 
to a single antiferromagnetic chain with next-nearest- 
neighbor interactions. In this paper we consider the spin- 
1/2 antiferromagnetic zigzag ladder with .anisotropy. The 
isotropic case has been studied beforeElEJ as a function 
of the coupling parameter, j = J2/J1, which is the ra- 
tio of the next-nearest-neighbor interaction, J2, to the 
nearest- neighbor interaction, J\. As j increases, the sys- 
tem goes from gapless (single chain) to a dimer phase and 
then to a spiral phase, where the structure factor has a 
maximum at a momentum 7r/2 < q < ir. The system 
has a spin gap in these last two phases, and it there- 
fore only displays short-range order. In the limit that 
the intra-chain interaction is much larger than the inter- 
chain interaction (J — ► 00) the two chains decouple and a 
gapless single chain behavior is recovered. It has been ar- 
gued that this only happens, strictly speaking, at j = 00: 
the spin gap becomes exppaentially small as j grows, but 
it remains non- vanishing .113 Recently, on the other hand, 
it has been proposed that incommensurate quasi-long- 
range spin correlations should be observed if easjtplane 
(XY) anisotropy is included in the zigzag ladder £3 This 
is argued to be due to the presence of a "twist" term 
that results from the inter-chain interaction. It has been 
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proposed that there is one gapless mode and one mode 
with a gap in the regime of strong XY anisotropy in the 
inter-chain coupling. Another prediction of this work is 
the existence of spontaneous local spin currents. This, 
however, has been refuted in ref. 16. Also, other recent 
numerical workli3 has failed to confirm the gapless nature 
of the groundstate in the anisotropic XY case at weak 
interchain coupling. Recent Density Matrix Renormal- 
ization Group (DMRG) resultstS suggest, however, that 
the zig-zag ladder does indeed show a gapless chiral phase 
as predicted in ref. 15. 

In this paper we use_exact diagonalizations, the W>d- 
ified Lanczos methodic and the Davidson methodcJ to 
address the possibility of a transition from a spin-gap 
regime to a gapless regime as a function of j when 
anisotropy is present in small 5 = 1/2 antiferromagnetic 
zig-zag ladders. We compute various probes to identify 
the different phases and study their behavior close to the 
transition points. We study in particular the spin stiff- 
ness, the dimer correlation function, the structure fac- 
tor, and analyze in detail the spectrum in the various 
parameter regimes. Since the zig-zag ladder effectively 
has both nearest-neighbor and second neighbor interac- 
tions, a stiffness tensor is required to account for these 
two types of interactions. The eigenvalues of this 2x2 
matrix then become the natural spin rigidities that we 
use to clarify the behavior of the system in the vari- 
ous regimes. The stiffness of a system is a particularly 
good measure of the long-range nature of the ground- 
state. Introducing twisted boundary conditions leads to 
a response in the energy if the quantum states are ex- 
tended (gapless case). On the other hand, the energy is 
insensitive to a change in the boundary conditions if the 
quantum states are localized (spin-gap case). Therefore, 
the stiffness with respect to such a twist is positive if the 
system is gapless and it is zero if the system has a gap. 
Also, the dimer correlation function naturally signals the 
dimer phase while the structure factor is a natural way 
to detect and study the spiral phase. 

Our results are consistent with a gapless excitation 
spectrum in the case of XY anisotropy at weak inter- 
chain coupling. We obtain qualitatively similar results 
for the isotropic case, however. This is most likely a fi- 
nite size effect due to the exponentially small spin gap 
that persists in the isotropic zig-zag at weak coupling in 
the thermodynamic limit S3 In fact, we show that the 
phase diagram that is obtained from the analysis of the 
spectrum for finite systems may be consistent with the 
field theory predictionEij after performing extrapolations 
to the thermodynamic limit. 

The paper is organized as follows: In section II we 
present the model and the quantities to be calculated. 
In section III we present our results and in section IV we 
summarize the work. Technical details concerning ex- 
act diagonalization are given in the Appendix. We also 
briefly review the extrapolation technique to the thermo- 
dynamic limit here. 



II. MODEL AND PROBES 

The anisotropic zigzag ladder is defined by the Hamil- 
tonian 

H = \j? Y £ (S+Sr +1 + S r S + +1 ) + J{ £ SfS? +1 

i i 

+ \ J 2 Y E ( s t s i+2 + srs+ +2 ) + jf E s*S! +2 . (i) 

i i 

The spin operators refer to spin 5=1/2 states, while 
the summation i — 1, ...,N runs along the "rib" of the 
zig-zag ladder. We shall parameterize the interactions 
by the coupling parameter j — J* Y j J* Y and by the 
anisotropy parameter J\j J* Y = A = J|/J^ y . (The 
isotropic case reduces to j = J2/J1 and A = 1.) We 
will set J± = 1 henceforth. Consider first the nearest- 
neighbor Heisenberg chain with anisotropy, which corre- 
sponds to both the weak-coupling (Ji = 0) and to the 
strong-coupling ( J 2 = 0) limits of the zig-zag ladder (see 
Fig. 1). The spectrum is gapless for the case of XY-. 
anisotropy, |A| < 1, as shown by the Bethe ansatz.Eil 
The excitation spectrum consists of spin-1/2 particles 
dubbed spinons. Since flipping one spin represents a spin- 
1 excitation, the spinons can only be created in pairs. 
Therefore the conventional spin 1 magnons are decon- 
fined into spin-1/2 spinons that propagate incoherently. 
In the regime where A < — 1, the groundstate is ferro- 
magnetic. When — 1 < A < 1 the leading spin configura- 
tion is the Neel state with the staggered magnetization 
lying within the XY plane. At A = 1 the groundstate 
is again in a Neel state, but with a staggered magnetiza- 
tion that can point in any direction. Last, the spectrum 
shows a gap in the Ising regime at A > 1. The ground- 
state, on the other hand, displays strict long-range Neel 
order, with the staggered magnetization directed along 
the z axis. 

We shall begin our study of the antiferromagnetic zig- 
zag ladder by analyzing the classical limit of the isotropic 
Heisenberg case first: J* Y = J± = Jf and J% = J11 — 
Jf as S — ► 00. A spiral state Sf — Se l6i yields an energy 
per site of E{6) = S 2 J\\ cos(26>) + S 2 J ± cosO. This mag- 
netic energy is minimized at a pitch angle 0q that satis- 
fies cos(?o = — 1 J1/J11 for inter-chain exchange couplings 
that are below a critical value J^ = 4 J|| . A ferromagnetic 
state on each chain occurs, on the other hand, at strong 
coupling Jx > Jjl, with a pitch angle of 9q = n. The 
spins are thus arranged antiparallel in between chains. To 
summarize, the system is in a spiral phase for J± < 4 J|| , 
while it is in a collinear phase for J± > 4Jy. The same 
holds true when only XF-coupling exists. In the case of 
Ising coupling only, on the other hand, we have the ef- 
fective model H = J|| J2i SfSt+2 + J± SfSf+v There 
are two possible groundstates. The first is the collinear 
state defined by Sf = ^, for i even and Sf = —5 for 
i odd (this has a degeneracy 2) with an energy per site 
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of E = S 2 (Ju — Jj). The other state is the antiferro- 
magnetic one defined by ±5f = i(— l) 1 / 2 for i even and 
±Sf = g(— 1)( J±1 )/ 2 , for 2 odd (this has a degeneracy 
2 x 2), with an energy per site of E = — S 2 JS. We have 
an Ising antiferromagnet for J\ < 2JS and a collincar 
Ising ferromagnet for J*[ > 2Jif, with a first order tran- 
sition separating the two phases. 

Consider now the stiffness tensor in the classical limit. 
Imposing a spiral spin configuration on the zig-zag ladder 
with a pitch angle 9, the energy per site of the classical 
J\ — J 2 model is then given as above by 



e = — = S 2 Ji cos( 
N 



+ 9 1 ) + S 2 J 2 cos{29 + 8 2 ), (2) 



where we have added small twists 9\ and 9 2 to the 
nearest- neighbor and the next-nearest-neighbor terms, 
respectively. For 6\ = = 9 2 , we have that cos6*o = — I 
for J\ > 4J2 and that cos^o = —4(^1/^2) f° r <A < 4J2, 
as stated above. The spin currents are then given by 
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= -5 2 Jisin6»o 



= -S 2 J 2 sin 26» 



(3) 



and the rigidity components by 



P11 
P12 
P22 



dji 

oe 1 

dji 
89 2 

oe 2 



= -S 2 Ji cos 6»o 



-S 2 J 2 cos 26> 



(4) 



Note that both the spin currents and the stiffnesses are 
independent of the anisotropy parameter A in the clas- 
sical limit: the isotropic and the XY anisotropic cases 
give the same results. In the collinear phase at J% > AJ 2 , 
the spin currents vanish (j\ = = j 2 ) and pn = S 2 J\, 
P12 = and P22 = —S 2 J 2 . On the other hand, in 
the spiral phase at J\ < 4J 2 the local spin currents 
are non- vanishing: ji = ±S' 2 J 1 [1 — ( J1/4J2) 2 ] 1 / 2 and 
j 2 = However, the total spin current j s = j\ + 2j 2 

is null. The stiffness tensor of this spiral phase is given by 
p u = \S 2 {J 2 /J 2 ), Pl2 = and p 22 = S 2 [J 2 - l(Jf/J 2 )}- 
The natural stiffness associated with the total spin j s is 
the response to an external twist that satisfies 9 2 = 20±, 
and is given by p s — pn + ^p\ 2 + ip 22 . It reduces 
to p s = S 2 (Ji — 4J2) in the collinear phase and to 
p s = 5 2 [4J2 — j(J 2 /J2)] in the spiral phase. These re- 
sults are displayed in Fig. 2. Here it is shown that p s is 
always positive and that vanishes at the classical transi- 
tion point between the collinear ferromagnet to the spiral 
phase. Recall that the nature of the groundstate changes 
across this transition. Also, we remark that in the spi- 
ral phase only the stiffness p s and the total current j s 



are "well behaved". The other components show spon- 
taneous spin currents, while p 22 is not always positive 
(we will return to this point later while discussing the 
quantum case). 

In the general quantum case we calculate the stiff- 
ness in the standard wayE3 We consider the Hamiltonian 
(I) with periodic boundary conditions imposing uniform 
twists around the z axis: 

H(e 1 ,e 2 ) = \j? Y Y.( s i s ^ el9l + s * s t + i e - 101 ) 

i 



i+2- 



(5) 



Here 9\ and 9 2 are two independent twists that act sepa- 
rately along the interchain and intrachain directions, re- 
spectively. Expanding the exponentials to second order 
we obtain the form 

H(9 l7 9 2 ) = H (0,0) + 9^ Jl +^Y, J f 

i i 

-^ 2 E^-^ 2 E^ 2 (6) 



where 



J? = 2'^^ (^tSi+2 — ^i+2) 



(7) 



are the spin currents along the interchain and intrachain 
directions, respectively and where 



(8) 



are the kinetic energy operators. Using second-order per- 
turbation theory, we then obtain that the groundstate 
stiffness tensor, p a p — (d 2 E/d9 a d9p)\o, is given by 

i 

- 2 E^^E <0 ^>><Hjfio>, (9) 

vj&O V i,j 

where |0 > is the groundstate and \v > are the ex- 
cited states (a, (3 = 1,2). The groundstate is assumed 
to be non-degenerate. Both the spin current operator 
J a = J2i J? an d the kinetic operator T a = J^i T? com- 
mute with the translation operator, T, and conserve total 
spin S z (where a = 1,2). Therefore, the states \v > in 
the stiffness formula (9) are the excited states within the 
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subspace of a given magnetization and momentum that 
contains the groundstate. 

As in the classical case discussed earlier, the stiffness 
can also be calculated directlyo taking numerical deriva- 
tives of the energy with respect to small twists, &\ and 82, 
that are imposed on the system. This procedure requires 
care that 8\ and 82 are small enough so that there be no 
level crossings. (We prefer to use the correlation func- 
tion method, even though the results using both meth- 
ods agree very well.) The change in energy due to small 
twists 0\ and 8 2 takes the form 

SE = \pnQ\ + P12M2 + \piiQ\ (10) 

in the absence of spontaneous spin currents. It is then 
natural to consider the eigenvalues of the stiffness tensor 

p± = \ (pa + P22) ± (11) 

and the determinant D = p+p- = pwpn — P\2- These 
eigenvalues will be computed using expression (9) for the 
stiffness tensor in the next section. 

We shall also calculate the correlation function asso- 
ciated with the dimerization via linear response theory. 
Imposing a small explicit dimerization, we consider the 
Hamiltonian 

h{8iM = \j? y x; (1 + (-0**0 (s+sr+i + srs+ +1 ) 

i 

+ l J 2 Y E( 1 + H)^) (S+Sr +2 + S r S + +2 ) 

i 

i 

+ ^E( 1 + H) i fc)s?s? +2 . (12) 

i 

Once again using second order perturbation theory we 
obtain that the susceptibility Xap — — {d 2 E/d5 a d8p)\o is 
given by 

X^2Re£^XMr +J 

< 0|(lf + M?)\v >< v\(T? + M/)|0 > 

(13) 

where 

Ml = j;sfs; +1 

Mf = 4SfSt +2 . (14) 

In the dimer correlation function, the factor (— l) l+J im- 
plies that the states \v > are contained in the same mag- 
netization subspace, but in the k — fco ± 7r momentum 
subspace, where fco is the groundstate momentum. Last, 
the structure factor is defined in the usual way: 



S(q) = X e lqr C r (15) 

r 

where the correlation function C r is defined by 

C r = ^^j:<0\S,S i+r \0>. (16) 

and is normalized such that the local correlation function 
(r = 0) is unity. 



III. RESULTS 

We now proceed to study the XX Z zig-zag model (1) 
using exact diagonalization of finite systems with sizes 
up to N — 28. The full energy spectrum is obtained for 
the smaller system sizes, N < 16, while only the ground- 
state and the first few excited states can be determined 
for the larger system sizes. The eigenvectors and eigen- 
values of Hamiltonian eq. (1) are then substituted into 
Eqs. (9,13,15) to compute the various correlation func- 
tions. Our main aim is to study the transitions between 
the various phases. 

Spectrum. Let us first survey the energy spectrum that 
is displayed by these small systems. We shall keep track 
of important quantum numbers associated with each en- 
ergy level, such as the momentum along the rib of the zig- 
zag, the spin and the parity. We shall also identify points 
in parameter space where low- lying levels cross, and use 
this to identify phase transitions in the system. This pro- 
cedure is known to yield accurate transitk*. points when 
applied to even relatively small systems. BE] 

Let us begin by determining the quantum numbers 
of the groundstate as a function of the size N for the 
S = 1/2 zig-zag antiferromagnet. Periodic boundary 
conditions are imposed throughout. The groundstate 
is a spin singlet in general due to the antiferromag- 
netic interactions. For strong enough coupling between 
chains, j = J2/J1 < 1/2, it has either momentum tt for 
iV = 4n+2 or momentum for N = An. For weak enough 
coupling between chains j > 1/2, on the other hand, the 
momentum oscillates between and tt as a function of the 
coupling parameter j and of the system size NH There 
are several points along j in this regime where the cor- 
responding energy levels for these two momentum values 
cross. The groundstate is degenerate at these points, and 
this is reflected by peaks in the dimer correlation func- 
tion (see Fig. 8). Such level crossings grow in number 
as the system size grows, and this indicates that the two 
singlet states in question are in fact degenerate in the 
thermodynamic limit. By the Lieb-Schultz-Mattis theo- 
rempJ this is consistent with a spin gap in the excitation 
spectrum that survives the thermodynamic limit in the 
weak-coupling regime j > 1/2. 

Consider now the specific case of an XXZ zig-zag 
chain (1) with N — 16 sites under periodic boundary 
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conditions (see Table I). In the isotropic case, A = 1, 
the states are organized into spin multiplets due to the 
SU(2) spin invariance. Again, the antiferromagnetic in- 
teractions imply that the groundstate is a spin singlet 
(S z = 0) in general. The system has three well defined 
regimes: (a) strong-coupling, (b) intermediate coupling 
and (c) weak-coupling. In the strong-coupling limit (a), 
j = J2/J1 —* 0, the groundstate has momentum k — 0. 
The first excited state forms a spin triplet in such case, 
with k = 7r, while the second excited state is another spin 
singlet with S z — and k — ir. At j = j*[ ~ 0.24 there 
is a level crossing where the first excited states and the 
second excited state interchange (here is the value of 
j at which the level crossing occurs for the system with 
size N). In the thermodynamic limit the two lowest sin- 
glet states (S z = 0, k = and S z — 0, k — ir) become 
degenerate and there is a finite gap to the next excited 
state (the triplet state S z — 0,±1, k = ir). Although 
the system begins to dimerize at this stage (6), antiferro- 
magnetic correlations along the rib of the zift-zag remain 
dominant up to the Majumdar-Ghosh pointO at j = 1/2. 
The groundstate is doubly degenerate for any system size 
at this point, where the two states are perfectly dimer- 
ized. Antiferromagnetic correlations within chains of the 
zig-zag become dominant beyond this point at weaker 
couplings j > 1/2 (see Fig. 9). Another level crossing 
occurs as j increases to about ~ 1.6 such that the 
first excited states become two triplets with S z = 0,±1, 
k = ±7r/2. The groundstate displays quantum numbers 
S z = and k = at this stage (c), and it is no longer 
degenerate. This remains so as j —* 00. It should be 
mentioned, however, that White and Affleck predict that 
a non-zero spin gap persists in the thermodynamic limit 
at weak-coupling j > 1 between chaias, and that it be- 
comes exponentially small as j growsE3 (we will return to 
this point later). 

The spectrum of the anisotropic 5=1/2 XX Z zig-zag 
ladder has also been studied previously in the strong- 
coupling pegime up to the Majumdar-Ghosh line (0 < 
j < 1/2)0 A gapless regime occurs for XY anisotropy 
A < 1 and strong coupling j < j c i(A), an Ising antifer- 
romagnet along the rib of the zig-zag that shows a spin 
gap in the excitation spectrum occurs for A > 1 and 
j < j c i(A), and a dimer phase regime that also has a 
spin gap exists at j > j c i(A) and any A (here j c i(A) 
is the transition line obtained from extrapolation to the 
thermodynamic limit). The line j = j c ±(A) separates the 
gapless phase from the dimer phase for A < 1, while it 
separates the dimer phase from the (Ising) Neel phase for 
A > 1. The line at A = 1 and j < j c i separates the XY 
gapless phase from the Ising phase. 

Consider again the specific case of N = 16 sites in this 
instance, with anisotropy parameters A = 0orA = 0.5 
(see Table I). The groundstate is a singlet with S z = 
and k = as j increases from the strong coupling limit 
at j = up to j = 1/2. The first excited state has de- 
generacy 2, with spin S z = ±1 and momentum k = tt 
inside this regime. The next excited state is again a sin- 



glet with S z = and k = tt. As j increases there exists 
a line of points, j — j^(A), such that the energy level 
of the excited state with S z — ±1, k = tt crosses the 
energy level of the excited state S z = 0, k = ttE3 For 
3 > ici(A) the two lowest states are the two singlets 
with $2 = 0, k = and S z = 0, k = tt. Again, these 
two states become degenerate (with a gap to the first 
excited state) in the thermodynamic limit. The ground- 
state must be a spin singlet with S z = due to the an- 
tiferromagnetic couplings. This excludes the possibility 
that any S z = ±1 state be degenerate with it. As a re- 
sult, the level crossing between the S z = ±1 and singlet 
states can be used as an accurate criterion to determine 
the phase transition between the gapless and the spin- 
gap regimes£l If we increase j up to j = 1/2, then the 
system is exactly degenerate for all system sizes. This is 
a feature of the Majumdar-Ghosh point, which has a per- 
fectly dimerized groundstate. The behavior of the system 
does not change much for intermediate coupling j > 1/2 
beyond this point, with the exception that the momen- 
tum of the two lowest states interchanges between k = 
and k = tt as j grows. At a larger value of j — be- 
tween 1.2 — 1.6, however, a level crossing occurs between 
the S z — 0, k — or tt state (first excited state) and 
a state with S z = ±1, k = ir/2 (second excited state). 
Notice that the momentum of the first excited state is 
now k = 7r/2. This is to be expected in this regime since 
the two chains are weakly coupled, and the periodicity 
doubles. For j > 2, in particular, the first excited state 
is now four-fold degenerate (S z — ±1, k — ±tt/2), and 
we might expect to fall back into a gapless regime since 
the first excited state is not a spin singlet. This level 
structure is in fact the same as for j — > 00. 

In Fig. 3 we present the crossing points between the 
states S z = 0, k = ir and S z = ±1, k — ±ir/2 that sig- 
nal this transition as function of the number of sites for 
several values of the anisotropy parameter A. We have 
attempted to extrapolate the crossing points, j^, to the 
thermodynamic limit by employing^standard (BST) al- 
gorithm due to Bulirsch and StoerBfl The method is ex- 
plained in the Appendix. The results of this extrapola- 
tion procedure are shown in Fig. 4, where j C 2 is plotted 
as a function of A. The curve j C 2 separates a spin-gap 
(dimer) phase from a gapless phase at small interchain 
couplings. As expected the value of j C 2 grows near the 
isotropic point. (It should tend. to j = 00 at A = 1 
according to White and AffleckB) 

Physical Probes. Correlation functions can also be 
used to determine the nature of different thermody- 
namic groundstates. The spin rigidity, in particular, can 
discriminate between phases that do and do not have 
spin gaps. The stiffness of the nearest-neighbor spin- 
1/2 Heiscnbecg,chain has been calculated exactly via the 
Bethe ansatz.EHlln the thermodynamic limit this solution 
yields 

PU TT . TT 1 
yTF = — Sin -, r- (17) 

1 n \ n / 



5 



for anisotropics 



cos ■ 



J! 

jXY 



(18) 



where n is a positive integer. This yields p\\j J^ Y = 1/4 
in the isotropic case and pu/ J* Y = 1/tt in the XY^&kSe. 
In Fig. 5 we plot the exact Bethe ansatz stiffnesses for 
the single chain and compare it with the numerical diag- 
onalization results for N = 8, 12, 16, 20 sites. The spec- 
trum is gapless from the F-limit up to the isotropic 
Heisenberg case. Beyond this critical point the spectrum 
acquires a gap due to the Ising anisotropy and the stiff- 
ness goes to zero. The stiffness remains positive beyond 
the transition point for the small system sizes that we 
studied, however. The anisotropy at which pn = J± /4 
extrapolates nicely to the Heisenberg point, A = 1, in 
the thermodynamic limit, N — > oo, however. 

The zigzag ladder (or the chain with next-nearest 
neighbor interactions), on the other hand, is not solvable 
by the Bethe ansatz, and so other methods become nec- 
essary to compute the spin rigidity. With this purpose in 
mind, exact diagonalization calculations of the isotropic 
spin- 1/2 Heisenberg chain with next-nearest neighbors 
interactions where carried out first by Bonca et al.E2l on 
small systems, where the diagonal component pn of the 
stiffness tensor was computed (see also ref. 31). We have 
completed this study by calculating the remaining com- 
ponents of the stiffness tensor (9), including cases with 
anisotropy. In Fig. 6 we show the results for the various 
stiffnesses for the zigzag ladder in the isotropic case as a 
function of j for N = 16, while the same set of results 
are shown for the XY-case (A = 0) in Fig. 7. In Figs. 
6a and 7a we show the stiffnesses pn, P22 and p±2- For 
small j we are in the limit of strong interchain coupling 
(small next-nearest neighbor interaction) and we recover 
the previous results:!^ p n > and P22 < 0. This is 
also found in the classical case. In the opposite weak- 
coupling limit of very high j, pn remains negative, but 
P22 becomes positive (as in the classical case). The lat- 
ter is consistent with the extreme limit of two decoupled 
antifcrromagnetic chains (Ji = 0). Also, pi2 becomes 
non-zero in the intermediate region where quantum fluc- 
tuations are stronger and where the transition between 
weak-coupling and strong-coupling occurs. Note that the 
classical analysis reveals that the fact that P12 becomes 
nonzero and the fact that pn becomes negative are purely 
quantum effects. 

In the initial exact diagonalization study for the 
isotropic case, Bonca et al. chose to identify the point 
at which the stiffness component p\\ turns negative with 
a transition point to a quantum disordered phase with a 
spin gap.EH Extrapolating to the thermodynamic limit 
they estimated the transition point to occur at j c = 
0.270 ± 0.005. We believe that it is a better idea to look 
at the eigenvalues (11) of the full stiffness tensor to de- 
termine possible phase transitions .0 In Figs. 6b and 7b 
we plot p + for different system sizes. This stiffness is al- 
ways positive for all j, while p_ is always negative. For 



small j, p+ is positive and of order unity. As j grows 
p+ becomes close to zero at a value of j that is close to 
the value where pn crosses zero, jf^. As j increases fur- 
ther p + again begins to grow appreciable near the point 
at which P22 turns positive, j^ 2 - A finite size analysis re- 
veals that as the size N grows the first "zero" of p+ (near 
j = j^i ) occurs at smaller values of j and extrapolates to 
a value close to the dimer transition j c \ determined by 
the level-crossing method. Similarly, we expect that the 
second transition to a nonzero stiffness p+ (near j = j^ 2 ) 
signals the transition to a gapless regime that extends 
up to j — 00, and that it might therefore also signal a 
decoupling transition. This happens, however, both for 
A = and for A = 1. We believe that this is due to a 
finite-size effect in the latter case as discussed above. The 
fact that finite-size effects in the spin stiffness of a single 
chain become larger as the XY anisotropy decreases (see 
Fig. 5) supports this claim. 

The level crossing at weak- coupling, j^, does not, how- 
ever, correlate well with jj^. The crossings defined by 
jp2 appear at smaller values of j as compared with j^ 2 - 
Although both have the same trend, apparently finite 
size effects are stronger in the calculation of the stiff- 
ness than in the determination of the level crossings. We 
have therefore limited ourselves to the extrapolation of 
the level crossing points j^. 

In Fig. 8 we show the behavior of the dimer corre- 
lation function \y\. For j > j£[ the dimer correlation 
function increases signaling the spontaneous dimerization 
observed in the thermodynamic limit. The various peaks 
signal level crossings of the groundstate between near- 
degenerate states with momentum k = and k = n as j 
varies. The Majumdar-Ghosh point at j = 1/2 is special 
since the groundstate is doubly degenerate for all system 
sizes in such case. Beyond j ~ 1 the dimer correlation 
function \n becomes small. (The susceptibilities X12 and 
X22 are always quite small). 

Finally, remnants of the spiral phase that exists in the 
classical limit, S — > 00, are clearly apparent in the mo- 
mentum dependence of the structure factor. In Fig. 9 
we show the structure factor as a function of momentum 
for several values of j in the Heisenberg case at A = 1. 
We see that the location of the maximum shifts from 
k m ax = 7T (for j < 0.5) to a value ir/2 < k max < ir when 
j > 0.5, thereby signaling the spiral phase. The results 
are similar for A = 0. 



IV. CONCLUSIONS 



The S = 1/2 antifcrromagnetic zig-zag ladder is a 
difficult problem to solve due to the intrinsic frustra- 
tion and to the criticality displayed by both the strong 
and the weak-coupling limits. The weak-coupling limit 
is particularly difficult in the absence of anisotropy, in 
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which case, antiferromagnetic correlations diverge expo- 
nentially.E3 This renders any numerical study of finite 
systems hard. 

In this paper, we have performed an analysis of the ex- 
act properties of such finite systems, looking at various 
correlation functions and the structure of the spectrum 
both in the isotropic and the anisotropic cases. We have 
looked at the spectrum and have computed the spin stiff- 
ness of the zig-zag ladder, and have thereby found evi- 
dence for a gapless regime at weak coupling that survives 
the thermodynamic limit in the case of XY anisotropy. 
However, the isotropic and the anisotropic cases look 
qualitatively similar. We believe that this is due to a 
strong finite-size effect in the former case. This claim 
is supported by the increase of finite-size effects in the 
stiffness of a single S — 1/2 antiferromagnetic chain with 
decreasing XY anisotropy, as shown in Fig. 5. 

It was previously showni that the dimer transition can 
be accurately determined in relatively small systems by 
studying the level crossing of the first and second ex- 
cited states after extrapolating to the thermodynamic 
limit. We have used a similar criterion to detect a possi- 
ble transition from the dimer phase to a gapless phase at 
weak-inter chain coupling. Using standard extrapolation 
techniques (see Appendix) we have constructed a phase 
diagram (see Fig. 4) in this regime that is in agreement 
with a recent proposal for ,ar|gapless spiral phase in the 
presence of XY anisotropyJla 

Also, we expect the spin stiffness p + to be a good 
measure of the nature of the spectrum. In particular, 
it can be used as an order parameter to distinguish gap- 
less from spin-gap phases. A positive stiffness indicates 
a gapless excitation spectrum and a null stiffness indi- 
cates a net spin-gap. The antiferromagnetic zig-zag for 
spin S — 1/2 showed an appreciable positive stiffness in 
the limit of strong interchain coupling (similar to the sin- 
gle chain case), a very small yet positive stiffness in the 
intermediate (spin-gap) regime, and then again an ap- 
preciable positive stiffness in the limit of weak interchain 
coupling (similar to two decoupled chains). This was true 
for all values of the (XX) anisotropy, A < 1. It is known 
from previous workJ13 however, that a spin gap is al- 
ways expected to be present in the isotropic case at weak 
coupling. The stiffness should therefore remain zero in 
the thermodynamic limit at A = 1 in the weak-coupling 
regime. We believe that the discrepancy between this 
expectation and our results is a strong finite size effect 
in the isotropic case that is due to thepexponentially di- 
verging antiferromagnetic correlations .E3 Clearly, larger 
systems need to be studied. 
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APPENDIX: 

Exact Diagonalization. The size of the Hilbert space 
under consideration can be considerably reduced using 
the symmetries of the problem. The Hamiltonian com- 
mutes with the total spin operator S T , with the transla- 
tion operator T, the spin flip operator R, and the space 
reflection operator L (i —> N+l — i,i = 1,..., N). In the 

absence of anisotropics the Hamiltonian also commutes 

- 2 

with the total spin operator (St) and the energy levels 
come in spin multiplets. 

The operator R commutes with T but, although L 
commutes with R, it does not commute with T, since 
LT = T- 1 L. 

The action of the local operators Sf and Sf is sim- 
ply given in the direct product basis |mf > ■ • • \m z N >, 
which are eigenvectors of S T . In general, these states 
are not eigenvectors of the additional symmetries. We 
consider first the translations, obtaining the classes of 
states which are closed under them. One starts with 
a state \a > and applies N a times the translation op- 
erator T until one finds T Na \a >= \a >, where N a is 
necessarily a divisor of N. The state \a > is the rep- 
resentative of this class. In combinatorial theory this 
is called a necklace and its periodic part, of length N a , 
a Lyndon word. The other classes are formed proceed- 
ing in the same manner starting with other states, not 
already used, until all the states have been exhausted. 
In each class, since T Na — 1, the possible eigenvectors 
of T are tN a (k) = e N <> , with k — 0, . . . , N a — 1, cor- 
responding to the momentum pk = The opera- 
tors P N (k) = Ei=o'~ 1 (t^fey) i are projectors, i.e, they 
satisfy PN(k)PN(k') = Sk,k> PN(k). They also have the 
property TP^(k) = tjv(k)PN(k). The projector PN(k) 
acting on the states of a class with N a elements reduces 
to -^-P/v a (fc), if tjy(k) is one of the t^ a (k) and gives 
zero otherwise. The normalized eigenvector of momen- 
tum pk formed from the class of the state \a > is given 
by \k,a >= y / N a P]y a (k)\a >. Since the Hamiltonian 
commutes with S T and T, the subspaces of fixed magne- 
tization and momentum are invariant subspaces and it is 
important to consider each of these subspaces separately. 
A general state is written as a linear combination of the 
states \k,a > corresponding to that magnetization. The 
state \k,a > is represented by \a >. This allows us to 
reduce the size of the basis to the number of representa- 
tives. 

The spin reflection changes the sign of the magnetiza- 
tion. Therefore it is relevant only for the classification of 
the states when the magnetization is zero. IF R\a > is 
not in the class of \a > both eigenvectors with r — 1 and 
r = — 1 can be formed. IF R\a > is in the class of \a >, 
it must be obtained from this state by ^ translations 
(implying that N a must be even), since for half- integer 
spin it is linearly independent of \a >. The state \k,a > 
is an eigenvalue of R with r = 1 for k even and r = — 1 
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for k odd. 

The space reflection operator reverses the momentum, 
since one has LPjv(fc) = -Pjv(— fc)L. The state L\a > can 
be in a different class or in the same class. The states 
\k,a > and | — fc,a > are linearly independent, except 
for k — 0, with eigenvalue I = 1, or k = for N a 
even, with eigenvalue I = — 1. As a result, one sees that, 
for zero magnetization and t = ±1, it is possible to use 
all the operators T, R and L in the construction of the 
states. This allows a further reduction of the size of the 
subspace under consideration. Since some classes have a 
definite eigenvalue of R or L they are simply excluded 
if their values are not those of the state which we are 
studying. 

Since the complete Hamiltonian commutes with the 
operator T only transitions to states with the same 
momentum are allowed, even if separate terms in the 
Hamiltonian allow them. If H\a >= Y, b aTv(a ' b) \ b >> 
where p(a,b) is an integer and \b > is the represen- 
tative of a class, one finally finds < k, b\H\k,a >= T^ m ^ = T^™^ + 

^'b a \p^t ^ Ar «(^)) P< ' a ' & ^ wnere the summation is re- 
stricted to the allowed transitions^. 

The modified Lanczos methodta is very useful to ob- 
tain the ground state of an Hamiltonian. Restricting 
our attention to a subspace of fixed magnetization and 
momentum we will find the ground state of the block 
Hamiltonian in that subspace. To obtain an approxi- 
mate ground state we choose a trial state 1-00 > that can 
not be orthogonal to the true ground-state. We define a 
state l^i > aslla 



Since the ground state and the first excited of the sys- 
tem have different momenta they are both ground states 
of different subspaces, and we apply the modified Lanc- 
zos method to these subspaces. If the system does not 
have anisotropies one can alternatively look for the com- 
ponent of the spin multiplet of the first excited state with 
magnetization 1. 

Extrapolation. The results for the infinite system can 
be estimated using-standard extrapolation methods like 
the BST methodo In the BST algorithm we look for 
the limit of a sequence of the type T(hj) = T + a,\h^ + 
a 2 /ij" + j = 0,- • • ,N P — 1 (N p being the number of 
data points), where hj is a sequence converging to zero 
as j — ► oo, corresponding to different system sizes Nj. 
Typically hj = with Nj = a + bj , for some values of 

a, b. The value of the m th iteration for the sequence is 
obtained from 



T, 



(j'm + l) 



T, 



0m) 



( si ) 



I - 



TTl-1 



- I 



(A7) 



with j m = 0, • • • , M(m) and where M(m) is the number 

of values of Tm at each iteration. It decreases by one, at 
each iteration, from M(0) = N p to M(N p — 1) = 1, when 
the iteration process is fulfilled. As initial values one de- 



fines T^l = and T W = T{hj). The extrapolated value 



-.00 



(i) 

Af„- 



T 



(0) 



101 >= 



g|0o > - < H > | 0o > 
(< H 2 > - < H > 



2^1/2 



(Al) 



is T^_ 1 and the estimated error is e = \T i 
Finally, uj is a free parameter which is adjusted such that 
the estimate of the error is a minimum. 



where < 0i|0i >= 1, < 0i|0o >= and < H n >=< 
0o | H n 1 0o >• Defining a matrix of the Hapiltonian in 
the basis |0o >, |0i > we can diagonalize itO obtaining 
a next order approximation for the energy eo, and ground 
state |0o >, with 



100 > = 



e = < H > +ba 

|0O > -H01 > 



(1 + a 



where 



and 



/ = 



b=(<H 2 >-<H > 2 ) 1/2 
a = f-(f + l) 1/2 

< H 3 > -3 < H >< H 2 > +2 < H > ? 



2(< H 2 > - < H > 



(A2) 

(A3) 

(A4) 
(A5) 

(A6) 



Taking |0 O > as the new |0 O > we can iterate the method 
to obtain a better estimate for the energy and the ground 
state. 
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Table I- Lowest energy levels for TV = 16. The states are represented by their S z values and momenta (S z ;k). 
When more than one state is represented this means they are degenerate. When the momenta are not 0, 7r or n/2 
the momentum of the state is represented by an integer, n, such that k — jj-n. 
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Figure Captions 



Fig. 1- Diagram of the zig-zag ladder. The nearest-neighbor interaction is parameterized by J\ and the next-nearest- 
neighbor interaction is parameterized by J2- 

Fig. 2- Classical stiffnesses pn/(JiS 2 ), P22HJ1S 2 ) and p s /(JiS 2 ) as a function of j = J2/J1. 

Fig. 3- Crossing points of the first excited state as a function of N for different values of A at weak coupling. 

Fig. 4- Phase diagram in the regime of weak interchain coupling. The critical coupling parameter j C 2 is obtained from 
the crossing of the first excited state with the second excited state after extrapolating to the thermodynamic limit 
(see the Appendix). 

Fig. 5- Comparison of the exact Bethe ansatz result for the spin stiffness of a single chain to the numerical results 
obtained from diagonalizing small systems as a function of A. 

Fig. 6- (a) Stiffnesses pn, pu and P22 as a function of j for A = 1 and N = 16. (b) p + as a function of j for various 
system sizes for A = 1. 

Fig. 7- (a) Stiffnesses pu, p\2 and P22 as a function of j for A = and N = 16. (b) p+ as a function of j for various 
system sizes for A = 0. 

Fig. 8- Dimer correlation function for N = 16 and A = 1 and A = 0. 

Fig. 9- Structure factor, S(q), for N = 16 and A = 1 for various coupling strengths. The peak is located at q = it 
in the single-chain limit, j = but this momentum tends to q = ir/2 as j increases. This indicates the presence of 
important spiral correlations in the weakly-coupled 5 = 1/2 zig-zag antiferromagnet. 
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S=1/2 Antiferromagnetic Chain 
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